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Narrow-gap semiconductor compounds like Pbi_ x Sn x Te and Hgi-xCd^Te present band 
inversion under compositional variation. In a band-inverted hetero junction the fundamental 
gap, defined as the difference between r 6 and T 8 energies, has opposite signs on each side 
|T||. Type III superlattices with band inversion of CdTe/HgTe and PbTe/SnTe has been 
successfully grown in the past [@,§J. One of the most conspicuous characteristic of band- 
inverted heterojunctions is the existence of interface states lying within the fundamental 
gap, provided that the two gaps overlap In IV- VI compounds those interface states 

are properly decribed by means of a two-band model using the effective k • p approximation. 
On the contrary, the analysis is more complex in II- VI compounds due to mixing with 
heavy-hole states since non centro-symmetry effects are not negligible in this case. The 
equation governing conduction- and valence-band envelope-functions in a simple two-band 
model, neglecting far-band corrections, is a Dirac-like equation. Exact solutions can then be 
found in view of this analogy because one can use elaborated techniques like those related to 
super symmetric quantum mechanics [J]. The aim of this paper is to present an alternative 
way of solution based on the so called point interaction potentials [P|JT0[] (any arbitrary sharply 
peaked potential approaching the 5-function limit) along with a Green function method. We 
believe that our treatment gives a very intuitive explanation of the origin of interfaces states 
while other approaches may obscure the way how those states arise. Moreover, the effects 
of external electric and magnetic fields can be included in a straightforward fashion, as we 
will show later. 

In the effective-mass approximation the electronic wave function is a sum of products 
of Bloch functions at the band-edge with slowly varying envelope-functions. The two-band 
model Hamiltonian in the absence of external fields is of the form 

H = v ± a y p y + v z a z p z + - (3E G (z), (1) 

where the Z axis is perpendicular to the heterojunction, Eg{z) stands for the position de- 
pendent gap, aty, a z and f3 are the usual 4x4 Dirac matrices, v± and v z are interband 
matrix elements having dimensions of velocity. As usual, it is assumed that these matrix 
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elements are constant through the whole heterostructure due to the similarity of the zone 
centre in both semiconductors. Since the gap depends only upon z, the transversal momen- 
tum is a constant of motion and we can set the Y axis parallel to this component. In the 
two-band case there are four envelope-functions including spin and we arrange them in a 
four component vector F(r). This vector satisfies the equation 



where V(z) gives the position of the gap centre. It is understood that the growth direction is 
[111]. The way V(z) changes from one layer to another is not well understood but, assuming 
that the interface states spread over distances much larger than the interface region, we can 
confidently consider it as a step-like function. Accordingly we take 



9 being the Heaviside step function. Here, the subscripts L and R mean left and right sides 
of the heterojunction, respectively. 

As we have already mentioned above, the momentum perpendicular to the interface is 
conserved, and therefore we look for solutions of the form 



HF(r) = [E-V(z)]F(t), 



(2) 



E G (z) = E GL 9(-z) + E GR 9(z), 



(3a) 



V(z) = V L 9(-z) + V R 9(z), 



(3b) 




(4) 



to Eq. (|2|). The function F(z) satisfies the following equation 





A simple way to solve this equation is the Feynman-Gell-Mann ansatz \) I 





After a little algebra we obtain 



( 



d 2 1 M 



dz 2 h 2 v 2 U 



E G {z) 2 -[E- V{z)] 2 + vlpl - i Aa z ((3 - X)5(z) = 0. (7) 
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For brevity we have defined 

A = T*- < 8a > 

A = 2 * J; , (8b) 

and we have used the relationships d6(±z)/dz = ±5(z). Note that in the case of band- 
inverted heterojunction E GR E GL < 0. 

It is worth mentioning that Eq. ([?]) is nothing but a Klein-Gordon equation with scalar- 
like and electrostatic-like terms depending on position (like a relativistic spinless particle 
with a position-dependent mass in an electric field as occurs in QED) plus a point interaction 
potential arising from the discontinuity of the gap and the gap centre. The occurrence of 
this short-range potential makes it possible the existence of bound states deep in the gap. 
In order to find the bound states we use a Green function formalism, similar to previously 
used in the case of the Dirac equation with point interaction potentials |12| . To this end, let 
us consider the Green function associated to Eq. ([7|) without the point interaction potential 



d 2 1 
+ 



\E G {zf-[E-V{z)f+vlpl 



G(z,z';E)=h5(z-z'), (9) 



\ dz 2 h 2 v 2 

where I4 stands for the 4x4 unity matrix. Since we are interested in bound state levels, 
boundary conditions reads 

Jim G(z,z;E)=0. (10) 

The Green function is a 4 x 4 matrix which permits the factorization G(z, z'\ E) = 
g(z, z'; E)I±, where g(z, z'; E) is a scalar function since the operator in the left hand side of 
Eq. (HD is scalar. 

The solution of Eq. (0) is then simply written out as follows 

X (z) = zAa z ((3 - A)Ho dz' G(z, z'; E)5(z') X (z'), 

= zAa z ((3-X)g(z,0;E) X (0). (11) 

It is assumed that the envelope-functions are continuous at the heterojunction so that the 
value x(0) is defined without ambiguities; this is completely different to what found in the 



Dirac equation for point intereaction potentials (see Ref. and references therein). Once 
the Green function is known, the 4- vector xi z ) can De obtain and using Eqs. (Q) and (|6]) the 
envelope-functions are finally determined. Bound state levels can be computed taking the 
limit z — > in Eq. ((TTJ) . To obtain nontrivial solutions we require the 4x4 determinant to 
vanish. Thus, using the definitions of A and A given in (B) we obtain 



1 



U 2 vl 



- A (E GR - E GL f - {V R - V L f 



1 



J7(0,0;£), 



(12) 



At this point we would like to stress that we just require the value of the Green function 
at the origin of the (X, X') plane if only the bound state levels are needed. In the absence 
of external fields, as we are considering here, this value is actually not difficult to obtain. 
Let u + and w_ be two independent, scalar solutions of the Sturm-Liouville problem (0) 
(dropping the point interaction term), vanishing at +oo and — oo, respectively. Therefore 
we can write 

g{0, 0; E) = — -, (13) 

W[u + ,u_\ 

where U-] is the Wronskian of the two solutions. Defining two real parameters 



Kr = ^ ^\eI r -(E-V r ? + vIp\, ( 14) 

the two independent solutions are u + = exp^—K^z) and U- = exp(KLz) so that g(0, 0; E) = 
2/(K R + K L ). Using Eq. (O) one finally obtains 



k r + k l = 1-J±(e G r - E GL y - (v R - v L y. (15) 

Kr and Kl should be real for obtaining exponentially decreasing envelope-functions as 
\z\ — ► 00 and then the gaps must overlap, i. e. (E GR — E GL ) 2 /4 > (Vr — Vl) 2 . This solution 
agrees with that previously proposed by Korenman and Drew ||. The reader is referred to 
Ref. for a fully discussion of its physical implications. Here we stress the main advantages 
of using our method. First of all, we have restricted ourselves to the case of no external 



potentials. Nevertheless, it is clear that applied electric or magnetic fields can be easily 
handled with minor modifications of the equations. Note that the crucial point is that one 
assumes that the Klein-Gordon equation without the point interaction potential arising from 
the abrupt interface can be solved exactly and the corresponding Green functions is explicitly 
written out. This is so for a large variety of electric and magnetic field configuration, as 
pointed out in Ref. (T^j. Thus, for instance, it is possible to investigate Landau levels 
in band-inverted hetero junctions in a rather simple way, instead of using more elaborated 



mathematical treatments, as those recently carried out by Aggasi |L4| . In addition, it is also 
possible to study confined Stark effect, a topic which remains open in the literature. The 
second aspect we remark is the fact that there is no need to use an abrupt heterojunction 
model, simulated by a step potential. The only requirement is that K^} and Kj} must 
be much larger than the interface itself, an implicit assumption when using the envelope- 
function formalism. Qualitatively the profile of the heterojunction is soliton-like |7j and, as a 
consequence, it derivative is a sharply peaked function. Thus the integral equation ( [LTD can 
be solved by a limiting process, in analogous way to the Dirac equation for sharply peaked 
functions approaching the ^-function limit |T2|] . To conclude, we feel that the approach we 
developed holds in a large variety of cases of practical interest and it may help in a better 
understanding of interface states in band-inverted heterojunctions. 
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